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We study a quenched SU(2) lattice gauge theory in 4d in which the spatial gauge ensemble {f/,} 
is generated from a 3d gauge-Higgs model and the timelike link variables are "reconstructed" 
from the Higgs fields. The resulting ensemble is used to study quenched quark propagation with 
non-zero chemical potential jJ.. While it proves possible to alter the strength of the inter-quark 
interaction by changing the parameters of the dimensionally reduced model, there is no evidence 
for any region of parameter space where quarks exhibit deconfined behaviour or thermodynamic 
observables scale as if there were a Fermi surface. 
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1. Introduction 

Lattice QCD suffers from the so-called Sign Problem; the Lagrangian density with Nf quark 
flavors has the form qMq; the functional measure det^'M(;U) = det^'M*(-^) imphes that for 
/X / the action is complex, rendering Monte Carlo importance sampling impracticable. 

Can one at least perform lattice QCD simulations in the quenched Nf ^0 limit, i.e. study the 
propagation of valence quarks with /i / through a non-perturbative gluon background [[T]]? 

Of course, the information extracted from such an approach could be at best qualitative, since 
(unlike the case of T > 0) the gauge field ensemble {U^} can only respond to / via virtual 
quark loops. Anyway, such information might be valuable in furnishing a non-perturbative defi- 
nition of the Fermi surface, whose existence is assumed in most phenomenological treatments of 
dense matter. 

In the context of a random matrix theory, Stephanov showed that the quenched theory 
should be thought of as the Nf ^0 limit of a QCD-like theory with not just Nf flavors of quark ^7 G 3 
of the SU(3) gauge group but also with Nf flavors of conjugate quark q'^ G 3. As a consequence, 
qq'^ bound states appear in the spectrum, resulting in baryons degenerate with light mesons. For 
/i/r » 1 there is an onset transition from the vacuum to a ground state with quark number density 
nq > . For QCD it occurs at /lo ~ (mn is the nucleon mass) but if conjugate quarks are 

present /lo ~ 

We try to modify the gluon background in some way so that color confinement no longer holds: 
then the role of qq'^^ excitations may not be so important in determining the ground state in the quark 
sector. Our hope is that valence quark propagation in such a background may qualitatively resemble 
that of the deconfined regime of the phase diagram at /i /T » 1 . We start with the 3d configurations 
characteristic of the deconfined phase found niT > Tc, pL /T \ produced by the approach to hot 
gauge theory known as Dimensional Reduction (DR). The quenched action is the ^d SU(2) gauge 
- adjoint Higgs model obtained by DR from 4d SU(2) given by Eqn. (4) of Ref . : 

•^3^ = i3 £ ( 1 - ^tr^^,'; ) +2£tr(<p;,<p^) -2jc^tr(<p^i7^,i(p^+fi7j,0 + A^(2tr((p.,(p,0 - l)^, 

x,i> 7 ^ ^ X x,i X 

where (p = ^(pa'Za represents the adjoint Higgs field. In the DR approach, all non-static modes of the 
gauge theory are integrated out leaving a 3d gauge-Higgs model describing the non-perturbative 
behaviour of the remaining static modes. 

Our goal is to study quark propagation through a non-confining quenched gluon background 
with IJ. y^O. Since chemical potential couples to quarks via the timelike component of the current 
HYYoY, this is an inherently four-dimensional problem. In order to generate such a background 
we take a 3d configuration generated by the DR simulation, motivated by the fact that it describes 
deconfining physics, and "reconstruct" the gauge field in the timelike direction via the prescription: 

=cos(f^/^;;^)+/--^^|=sin(f^/^;^), (i.i) 

with g = Spatial link variables Ui are taken to be time independent and identical to their 3d 
counterparts. Henceforth we use 4d configurations {U^} generated as outlined above as input in 
quenched studies. 
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Since q and q fall in equivalent representations of the gauge group hadron multiplets contain 
both qq mesons and qq, qq baryons, which are degenerate at /i = 0. In the chiral limit the lightest 
hadrons are Goldstone bosons and can be analysed using chiral perturbation theory (;^PT) [Q] : at 
leading order for » T a second order onset transition from vacuum to matter consisting of tightly 
bound diquark scalar bosons is predicted at exactly /lo = niji/l. In the limit pL }JLo+ the matter 
in the ground state becomes arbitrarily dilute, weakly-interacting and non-relativistic: a textbook 
example of Bose-Einstein condensation. At the same point the chiral condensate {qq) starts to 
fall below its vacuum value and a non-vanishing diquark condensate {qq) develops. The diquark 
condensate spontaneously breaks U(l) baryon number symmetry, so the ground state is superfluid. 

More recent simulations have found evidence for a second transition at larger /i to a decon- 
fined phase, as evidenced by a non-vanishing Polyakov loop |j5|] and by a fall in the topological 
susceptibility [^. In this regime thermodynamic quantities scale according to the expectations of 
free field theory (also referred to as "Stefan-Boltzmann" (SB) scaling), namely Uq oc jj.^, and energy 
density £ <=< 

2. Numerical Results 

We have chosen jS = 9.0 sufficiently close to the continuum limit for the DR formalism to 
be trustworthy, and start with a point with T = IT^: this corresponds to K" = 0.3620027, A = 
0.0020531. It is important to make a precise determination of the pion mass at = 0: we obtain 
rrijcat = 0.2321(1) on a 8^ x 32. Since this scale is not too dissimilar to Lj^, we have repeated 
the measurement on 16^ x 64, where we find mj^at = 0.2368(3). The systematic error due to finite 
volume is significant, but small enough at 2% to be acceptable for this exploratory study. We 
have analysed the quark density and chiral condensate {qq) as functions of jJ. for various j with 
< n/T < S. There is a transition at jxaf 0.12, becoming more abrupt as diquark source j — > 0, 
from a phase with riq = 0, {qq) constant to one in which Ug increases approximately linearly with 
/I and {qq) oc jj^j^ jj^ complete accordance with the scenario described by ;^PT in which 
as n increases at T « there is a transition at /i^ = mjt/2 from the vacuum to a weakly-interacting 
Bose gas formed from scalar diquarks. The diquarks are supposed to Bose-condense to form a 
superfluid condensate; on a finite system this must be checked at j 7^ using the {qq+) observable: 
the condensate increasing monotonically with jx. To determine the nature of the ground state an 
extrapolation 7 ^ is needed. We have used a cubic polynomial for data with 0.02 < 7a < 0.1; 
Fig. [l| confirms that once again there is an abrupt change of behaviour in the order parameter at /i 
mji/2, and that the high-/i phase is superfluid. Next we explored a parameter set corresponding to a 
smaller scalar "stiffness" by changing to K" = 0.1. Naively this corresponds to a huge value of T /Tc, 
ie. taking us further into the deconfined phase of the hot gauge theory. Of course, whether DR- 
based concepts remain valid for the reconstructed theory must be addressed empirically. This time 
we used a volume 8^ x 64 for the bulk observables and at = determined the pion mass nijiat = 
0.1377(1) on 8^ x 128, and nijiat = 0.1423(4) on 16^ x 64, showing that the finite volume error is 
now roughly 3%. We studied nq and {qq) as functions of /X for < /i/r w 12. It is noteworthy that 
for n > mji/2 riq^j.) is numerically very similar to the values found at K" = 0.3620027, whereas 
for n < niji/l the chiral condensate {qq) is significantly smaller, indicative of a weaker quark - 
anti-quark binding at this smaller K. Moreover, using the Gell-Mann-Oakes-Renner relation, valid 
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Figure 1: {qq+) vs. pLat for various j; the extrapolation to j ^ is obtained using a cubic polynomial. 



for xVY, flml = Im^iqq) we have: f^a^ = 0.99 for K = 0.1 and f^a^ = 0.58 for K = 0.3620027. 
The non-interacting limit is /;t — > oo therefore the interparticle intractions are weaker at K" = 0.1. As 
before, there is a clear discontinuity in the observables' behaviour at Hc — mjc/2, and the general 
picture is qualitatively very similar, suggesting that the ;tPT scenario is still applicable. Diquark 
binding is now also much weaker. 

It is disappointing that we have found no qualitative change in physics as the parameters are 
varied - recall that the ;^PT model which describes the results reasonably well is based on the 
assumption of confinement, or at least on the presence of very tightly bound diquark states in the 
spectrum. 

To explore the parameter space more widely we focussed on a single observable, n^, and 
scanned the plane on 8^ x 16 at five different values of /i with jS = 9.0, ma = 0.05 and 

ja = 0.01. Data shows that except for A = 0.1 the results for fixed n are practically independent 
of K and of A. Moreover, tiq increases linearly with jj. over a wide region of parameter space, 
as it does for jj. > Hc- This approximate linear behaviour is once again a prediction of ;tPT, and 
is to be contrasted with the riq <x jx^ behaviour expected of a deconfined theory where baryons 
can be identified with degenerate quark states occupying a Fermi sphere of radius kp ^ jx. The 
absence of this scaling is a further reason to conclude that the reconstructed model does not describe 
deconfined physics. 

2.1 Bosonic Spectrum 

An appropriate set of states to look at, to study the bosonic spectrum, is constituted of pion, 
scalar, higgs and goldstone. Pion and scalar states are related via the U(l)e global symmetry % ^ 
e"^^X, X ^ x^'"^- Analogous to chiral symmetry for continuum spinors, this is an exact symmetry 
of the action in the limit m ^ 0. In a phase with spontaneously broken chiral symmetry, the scalar 
is massive, and the pion a Goldstone mode, becoming massless as m ^ 0. Similarly, "higgs" and 
"goldstone" diquark states are related via the U(1)b baryon number rotation x ^ e'^^X' X ^ X,^ '^, 
an exact symmetry of the action in the limit j ^ 0. In a superfluid phase with {qq+) / 0, the higgs 
is massive, and the goldstone massless in the limit j —>■ 0. 

The boson correlators are constructed from the Gor'kov propagator where appear the 2 x 2 (in 
color space) matrices N ~ {XxXy) ^rid ^ ~ iXxXy) which are known as the normal and anomalous 
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Figure 2: Mass spectrum of various bosonic excitations as a function of pLa,. 



parts respectively. On a finite volume A = for 7 = 0; limj^olimy^ooA 7^ signals particle-hole 
mixing resulting from the breakdown of U(l)e symmetry, and hence superfluidity. Due to SU(2) 
symmetries the only independent components of are ReNw = N and ImAi2 = A and their barred 
counterparts. 

We have studied the model with a chemical potential jiat = 0.25, ie. above the critical pLc 
required to enter the superfluid phase. All four channels yield clear signals for single particle 
bound states. 

Like any meson constructed from staggered fermions, the correlators in principle describe two 
states and must be fitted using the form 

C{t) =A[e-"" + e-"'^^'-'^]+B[e-^' + {-\ye-^^^'-% (2.1) 

where m and M denote the masses of states with opposite parities. In most cases we find M ^ m; 
however for /i > the pion correlator has a distinct "saw-tooth" shape, and in fact the fit yields 
riijc > Mhi, where n denotes the usual pseudoscalar pion, and bl a state of opposite parity, which 
must therefore be scalar. In Fig. § the corresponding spectrum for all the boson states. Note that all 
states are approximately degenerate at /i = 0; the equality of pion, higgs and goldstone correlators 
is guaranteed by SU(2) symmetry at = 0, but the degeneracy of the scalar in the chirally-broken 
vacuum can only arise as a result of meson-diquark mixing due to j / 0. Next, note that the pion 
mass remains constant for jj. < jj.^, where it is a pseudo-Goldstone boson associated with chiral 
symmetry breaking, and then falls once the superfluid phase is entered; this is in accordance with 
the predictions of ;tPT. Most of the other states show a much steeper decrease with jj. for jj. < jj.^, 
followed by a gentle rise to a plateau at ma 0.13 in the superfluid phase jj. > /i^., precisely as 
expected of the goldstone state in the superfluid phase with diquark source j / 0. 

The exception is the higgs, which rises more steeply to become the heaviest state at large /i. We 
conclude that the breaking of degeneracy between higgs and goldstone states is clear supplementary 
evidence for the breaking of U(1)b symmetry in the superfluid phase and all states with = 0+ 
including the bl but except the higgs have some projection onto the Goldstone state, regardless of 
whether the original interpolating operator is mesonic or baryonic. 

2.2 Fermionic Spectrum 

We have also studied the fermion spectrum, often in the context of condensed matter called 



5 



Two Color Matter in the Quenched Approximation 



Pietro Giudice 



lattice 32x8x8x64 
am=0.05 
j=0.02 
k =0.0 



o 


n.=o.o 


□ 


n=o.i 




^=0.2 


A 


H=0.3 


< 


H=0.4 



(3=9.0 

K=0.3620027 
?i=0.0020531 



® ® ® 



10 20 30 40 50 

t 



Figure 3: Close-up of A(f)for various /xa,. 
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Figure 4: E and Y vs. ^at in both normal and 
anomalous channels. 



the quasiparticle spectrum. Since the Gor'kov propagator ^ is not gauge invariant we have to 
specify a gauge fixing procedure. A feature of the quenched approach is that it permits large 
statistics to be accumulated with relatively little CPU effort. This has enabled us for the first time 
in a gauge theory context to study W at jj. 0, by helping to overcome the sampling problems 
associated with gauge fixing. We have experimented with two gauge choices: Unitary gauge 
(p ^ (p' = (0,0,93), which is implemented before the reconstruction of the 4th dimension, and 
is unique up to a Z2 factor, specified by demanding 93 > 0; and Coulomb gauge, implemented 
by maximising Y.xi^^{^xj + ^x-n)' attempt to make the gauge fields as smooth as possible 
and hence improve the signal-to-noise ratio. We have studied the normal and anomalous fermion 
timeshce propagators on a 32 x 8^ x 64 lattice at j8 = 9.0, k: = 0.1, A = 0.0020531, ma = 0.05, 
ja = 0.0.2 and jxat = 0.3; the last value chosen to ensure /i > /if. Two features to note are that the 
Coulomb data is roughly twice as large as the unitary data reflecting an enhanced signal and that 
there is little variation with kx- 

In the NJL model the quasiparticle propagator can be successfully fitted using the forms 

N{t)= Pe^^"' + 2f^-'^'v(L,-0 ^ (2.2) 
A{t)=R[e-^^' -e''^^^^'-'^ (2.3) 

where for /i 7^ there is no reason to expect P = Q, but for a well-defined quasiparticle state the 
equality Ej^ = Ea should hold. Here, by contrast, only the anomalous channel fits produced an 
acceptable shows any evidence of gauge independence. The value of Ea obtained is very 

close to mji/2, indicating that at this value of K the pion is a weakly bound state. Another striking 
feature of the data is the approximate forwards-backwards symmetry ofN{t), implying P ^ Q. We 
have studied the dispersion relations SaI^.y) for data taken on a 32 X 8^ X 64 lattice; it confirms 
that the quasiparticle excitation energies are ^-independent: one motivation for our study was to 
investigate to what extent the concept of a Fermi surface, which is not strictly gauge invariant, can 
be put on a firm empirical footing in a gauge theory. Our results shows no evidence for a Fermi 
surface. The nature of the quasiparticle excitation is clarified a little at the other parameter set 
studied, namely K = 0.3620027. In this case our results show no evidence for any well-defined spin- 
i state in either normal or anomalous channels; as a result of confinement the excitation spectrum 
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of the model seems to be saturated by the tightly bound spin-0 states of Fig. 0. Fig. || shows a close- 
up ofA{t) for various pL values, showing the presence of an oscillatory component whose amplitude 
initially grows with /i (the pLat = 0.3 points overlay those from jxat = 0.2), but whose wavelength 
is roughly -independent. The origin of the oscillation could possibly be associated with the non- 
unitarity of the model, but is most likely a manifestation of independent spin-^ excitations being 
ill-defined due to confinement. 

We have fitted the k = 0.3620027 data to the forms 

N{t) = Pe-^''' co^{TNt + + Qe-'^"^^'-*'^ co^{TN{Lt - t) + <^), (2.4) 
A{t) =R[e-'^^'cos{rAt + (t))-e-^^^^'-'^cos(rAm-t) + (!))], (2.5) 

where we interpret E as the energy and F as the width of a quasiparticle excitation. The results' 
most striking feature is their independence of /i, with F of the same order of magnitude as E. 
An interesting systematic effect is that En > Ea while Fat < Fa, which has motivated us to study 
Vf'^ + F^ vs. }X : the disparity between normal and anomalous channels is significantly reduced. 
Inspection of the /ic?, = 0.3 data also shows that the gauge dependence of this result is 0(20%) 
at worst. Numerically, VE^ + > nin, indicating strong quark - anti-quark binding, due to the 
persistence of confinement at this value of K. Therefore we can interpret the effect of confinement 
as rotating the quasiparticle pole into the complex plane, the rotation angle being larger in the 
anomalous channel than in the normal one. 

3. Conclusions 

Our attempt to alter the nature of the gluon background by changing the parameters of the 3d 
DR gauge-Higgs model has been a partial success, in that in going from k = 0.3620027 to = 0.1 
the strength of the binding between quarks weakens significantly. 

However, in neither case is there evidence for significant departure of Uq, (qq) and (qq) from 
the behaviour predicted by ^PT, so that even if quarks are important degrees of freedom at K" = 0.1, 
there is no evidence for the formation of a degenerate system signalled by SB scaling. 

Sadly though, it appears to remain the case that despite its "unreasonable effectiveness" in 
virtually all other aspects of lattice QCD, the quenched approximation has nothing useful to tell us 
about the physics of high quark density. 
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